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, Abstract 

Radio networks (RN) are distributed systems (ad hoc networks) consisting in 
n > 2 radio stations. Assuming the number n unknown, two distinct models of 
RN without collision detection (no-CD) are addressed: the model with weak no-CD 
£ — . ' RN and the one with strong no- CD RN. We design and analyze two distributed 

leader election protocols, each one running in each of the above two (no-CD RN) 
models, respectively. Both randomized protocols are shown to elect a leader within 
C(log (n)) expected time, with no station being awake for more than C(loglog (n)) 
(-H ' time slots (such algorithms are said to be energy-efficient) . Therefore, a new class of 

c/3 . efficient algorithms is set up that match the r2(log (n)) time lower-bound established 



by Kushilevitz and Mansour in [11]. 



1 Introduction 

00 

^2 '• Electing a leader is a fundamental problem in distributed systems and it is studied in a 

t^- ' variety of contexts including radio networks [3]. A radio network (RN, for short) can be 

viewed as a distributed system of n radio stations with no central controller. The stations 
are bulk-produced, hand-held devices and are also assumed to be indistinguishable: no 
identification numbers (or IDs) are available. A large body of research has already 
focused on finding efficient solutions to elect one station among an n-stations RN under 
various assumptions (see e.g. [3, 11, 17]). It is also assumed that the stations run on 
batteries. Therefore, saving battery power is important, since recharging batteries may 
not be possible in standard working conditions. We are interested in designing power- 



saving protocols (also called energy- efficient protocols). The present work is motivated 
by various applications in emerging technologies: from wireless communications, cellular 
telephony, cellular data, etc., to simple hand-held multimedia services. 
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The models. As customary, the time is assumed to be slotted, the stations work syn- 
chronously and they have no IDs available. No a priori knowledge is assumed on the 
number n > 2 of stations involved in the RN: neither a (non-trivial) lower-bound nor 
an upper-bound on n. During a time slot, a station may be awake or sleeping. When 
sleeping, a station is unable to listen or to send a message. Awake stations communicate 
globally, i.e. the underlying graph is a clique, by using a unique radio frequency channel 
with no collision detection (no-CD for short) mechanism. In each time slot, the status 
of the unique channel can be in one of the following two states: 

• either SINGLE: there is exactly one transmitting station, 

• or NULL: there is either no station or more than two (> 2) broadcasting stations. 

When the status is NULL, each listening station hears some noise and can not decide 
whether or more than 2 station are broadcasting. When the status is SINGLE, each 
listening station hears clearly the message sent by the unique broadcasting station. 

In the weak no-CD model of RN, during a time slot each awake station may either 
send (broadcast) a message or listen to the channel, exclusively. By contrast, in the 
strong no-CD model of RN, both operations can be performed simultaneously by each 
awake station, during a time slot. Hence, in the strong no-CD model, when exactly 
one station sends at time slot t, then all the stations that listen at time t, transmitter 
included, eventually receive the message. In the literature, the no-CD RN usually means 
the strong model of RN, see e.g. [11, 14]. In the weak no-CD case, such a transmitting 
station is not aware of the channel status. 

Such models feature concrete situations; in particular, the lack of feedback mecha- 
nism experiences real-life applications (see e.g. [12]). Usually, the natural noise existing 
within radio channels makes it impossible to carry out a message collision detection. 

Related works. The model of RN considered is the broadcast network model (see 
e.g. [3]). In this setting, the results of Willard [17], Greenberg et al. [7] (with collision 
detection) or Kushilevitz and Mansour [11] (no-CD) for example, are among the most 
popular leader election protocols. 

In such a model, Massey and Mathys [12] serves as a global reference for basic 
conflict-resolution based protocols. Previous researches on RN with multiple-access 
channel mainly concern stations that are kept awake during all the running time of 
a protocol even when such stations are the "very first losers" of a coin nipping game 
algorithm [15]. In [9], the authors design an energy-efficient protocol (with o(loglog(n)) 
energy cost) that approximate n up to a constant factor. However, the running time 
achieved is O (log 2+e (n)) in strong no-CD RN. Also, an important issue in RN is to per- 
form distribution analysis of various randomized election protocols as derived in [4, 8] 
for example. 

Our results. Two leader election protocols are provided in the present paper. The 
first one (Algorithm 1) is designed for the strong no-CD model of RN, while the second 
one (Algorithm 2) is designed for the weak no-CD model of RN. Both are double-loop 
randomized algorithms, which use a simple coin-tossing procedure (rejection algorithms). 
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Our leader election protocols achieve a (quasi) optimal O(logn) average running time 
complexity, with no station in the RN being awake for more than (9(log log (n)) time 
slots. Indeed, both algorithms match the f2(logn) time lower-bound established in [11] 
and also allow the stations to keep sleeping most of the time. In other words, each 
algorithm greatly reduces the total awake time slots of the n stations: shrinking from 
the usual O(nlogn) down to O(nloglogn), while their expected time complexity still 
remains O(logn). Our protocols are thus "energy-efficient" and suitable for hand- held 
devices working with batteries. 

Besides the algorithms use a parameter a which works as a flexible regulator. By 
tuning the value of a the running time ratio of each protocol to its energy consumption 
may be adjusted : the running time and the awake duration are functions of a. 

Furthermore, the design of Algorithms 1 and 2 suggests that within both the weak 
and the strong no-CD RN, the average time complexity of each algorithm only differs 
from a slight constant factor. Also, our results improve on [13]. 

Outline of the paper. In Section 2, we present Algorithm 1 and Algorithm 2, and 
the main complexity results are given: Theorem 1 and Theorem 2, corresponding to 
Algorithm 1 and Algorithm 2. Section 3 and 4 are devoted to the analysis of both 
algorithms, by means of tight asymptotic techniques. 

2 Algorithms and Results 

Both algorithms rely on the intuitive evidence that all stations must be awake altogether 
within a sequence of predetermined time slots in order to be informed that the election 
succeeds of fails. 

To realize the election a first naive idea is to have stations using probabilities 1/2, 
1/4, ...to wake up and broadcast. This solution is not correct since with probability 
> never a station ever broadcasts alone. In order to correct the failure we will plan 
an unbounded sequence of rounds with predetermined length. Awake time slots are 
programmed at the end of each rounds in order to allow all stations to detect the possible 
termination of the session. 

In the following, we let a denote a real (tuning) parameter value, which is required to 
be > 1 (this is explained in Remarks 2 and 3). 

2.1 Algorithm 1 

In Algorithm 1 the stations work independently. Given a round j in the outer loop 
(repeat-until loop), during the execution time of the inner loop each station randomly 
chooses to sleep or to be awake : in this last case it listens and broadcast, simultaneously. 
If a unique station is broadcasting, this station knows the status of the radio channel; 
if the status is SINGLE, it becomes a candidate. At the end of round j, every station 
wakes up and listens to the channel and the candidates broadcast. If there is a single 
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candidate, the status is SINGLE again, and this candidate is elected. Every listening 
station knows the channel status, and is informed that the election is done. Otherwise, 
the next round begins. 

round <— 1; 
Repeat 

For k from 1 to [a roura(i ] do /* probabilistic phase */ 

Each station wakes up independently with probability l/2 k . 
An awake station listens and broadcasts 

If a unique station broadcasts Then it becomes a candidate station Endlf 
EndFor /* deterministic phase */ 

At the end of each round, all stations wake up, listen, 
and in addition all candidate stations broadcast; 
If there is a unique candidate then it is elected Endlf; 

round <— round + 1; 

until a station is elected 



Algorithm 1. Leader election protocol for strong no-CD RN 



The brackets in both algorithm represent the actions that take place in one time 
slot. Notice that the content of the broadcasting message is not specified, since it has 
no importance. The status "candidate" is valid for one round duration only. 

Definition and Notation 1 For the sake of simplicity, the following notations are 
used throughout the paper. We let 

f = j*(n) =' riog a log 2 (n)l. (1) 
Let also C :R 2 —>R defined by 

C(x,y) *T JEi- (2) 



(y-l)(l-y(l-x))' 



We have 



Theorem 1 Let a £ (1,1.1743...) and p\ = .14846... On the average, Algorithm 1 
elects a leader in at most C(p*, a) log 2 in) + O(loglogn) time slots, with no station 

being awake for more than 21og a log 2 (n) ^1 + o(l)j mean time slots. 
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2.2 Algorithm 2 



In the case of the weak no-CD model of RN a potential candidate alone cannot be aware 
of its status since it cannot broadcast and listen simultaneously. The awake stations 
choose to broadcasts or (exclusively) to listen. In the inner loop of Algorithm 2 any 
exclusively broadcasting station is called an initiator; when there is a unique initiator, a 
station that is listening to the channel, hears clearly the message. Such a station is said 
to be a witness. Witnesses are intended to acknowledge an initiator in the case when 
there is exactly one initiator. 

round <— 1; 
Repeat 

For k from 1 to [a round ] do /* probabilistic phase */ 

Each station wakes up independently with probability l/2 fc ; 
With probability 1/2, each awake station either broadcasts 
the message(/c) or listens, exclusively; 
If there is a unique initiator, the status is SINGLE; 
The witness(es) record(s) the value (k) Endlf 
EndFor /* deterministic phase */ 

At the end of each round, all stations wake up; 
All witnesses forward the recorded message; 
If the status is SINGLE 

there is a unique witness, the other stations receive the message (k) 
Then the unique initiator of the message (k) is elected 

[ and replies to all the stations to advise of its status Endlf 

round <— round + 1; 
until a station is elected. 



Algorithm 2. Leader election protocol for weak no-CD RN 



During the last deterministic phase all stations wake up. An election takes place in 
a round if this round is the first when there is a unique time slot in which there remain 
exactly a unique initiator and a unique witness. When a station is a witness, it continues 
to behave as it were not. Hence, a witness can be a witness twice or more, or even be 
an initiator. This is obviously not optimal but simpler for the analysis. 

We have 

Theorem 2 Let a G (1,1.0861...) and p\ = .07929... On the average, Algorithm 2 
elects a leader in at most C{j>\, a) log 2 (n) + O(loglogn) time slots, with no station 

being awake for more than 5/2 log a log 2 (n) ( 1 + o(l) J mean time slots. 
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3 Analysis of Algorithm 1 



Assume that Algorithm 1 begins with its variable round = j, and let pj be the probability 
that one station is elected at the end of this round j. One may also view pj has the 
conditional probability that an election occurs at round j knowing that it did not occur 
before. Within that round, that is for k ranging from 1 to \a J ~\ , every station decides to 
broadcast with the sequence of probabilities (l/2 h ) 1<k< ^ a jy 

The probability pu n ) that there is exactly one candidate for k = i in a round is 

/ i \ n—l 

n / 1 



p(i, n ) = 2? y- - 2? J • We then have 

\a>~\ \ai~\ \a*~\ 

pj = e^mIK 1- ^.")) = i 1 -^))" 1 ni 1- ^)) (3) 

k=l i =1 k=l i=l 

and pj = tj Sj , with 

l 3 = 12 (P(k,n)) (m+1) and Sj = IJ(l-P(i, n ))- ( 4 ) 

m=0 k=l i=l 

Note that Sj represents the probability of having no candidate in the jth round. 

Remark 1 Simple considerations show that when 2 aJ <C n, the probability Sj is close 
to 1 (and is "far from" 1 when 2 aJ 3> n). This remark explains the occurrences of the 
crucial values n/2 aJ and the definition of j* in Eq. (1). 

3.1 Two Lemmas for a Lower Bound on pj 

Lemma 1 provides a lower bound on the finite product denoted by Sj in Eq. (4). 
Lemma 1 Let j = j(n) such that j(n) — > +oo and n/2 Qj — > 0. We have 

liminfs,- > .1809.... 

n J 

Proof: We divide the product Sj in (4) as follows: for any r G {2, . . . , [a- 7 ]}, 

r- 1 \<* j ~\ 

s 3 = II^ 1 " P{i,n)) X II ( X " P(i,n))- ( 5 ) 
i=l i=r 

Since for any a € [0, 1] and n > 1, (1 — a) n < e~ ara , we have for any i > 1, 

1-P(i,„) > 1 - ^ exp ^— J > l--exp(--) e 1 ^; ( 6 ) 
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and more precisely, for i £ {r, . . . , [a-']}, with r large enough, 

> 1 -$ CT P(-5) exp (^) " 1 -? CTP (-?)( 1 + | : )- (7) 
Now we let r = r(n) = [1/2 log 2 (n)J , ensuring that 2 r /y / n is bounded. A simple 

r-1 

bounding argument in Eq. (6) shows that JJ(1 — P(i,n)) ~^ 1 when n — ► oo. 

i=i 

Thus, using also inequation (7), Eq. (5) yields 

limiof aj > exp In (l - (l + |) J e x p (-£)_ 

\ m>l i=l 

Denote by A n the right hand side of the inequality Computing the second sum (in- 
side the exponential in A n ) is completed through asymptotic approximations using the 
Mellin transform in Lemma 4 given in Appendix A (this is possible since n/2 a — > by 
assumption). By Lemma 4, 

A n > exp ( o(l) - £ (1 m + J^ ! - sup |tf„(log a n)| ] , (9) 

m>l 

where sup„ |[/^(z)| < .024234 for all z, as stated in Lemma 4 (Appendix A). 
Therefore, the third term in Eq. (9) is bounded from below by 



exp ( - .024234 ^ = .96092 . . . 



m>l 



Next, the second term within the exponential in Eq. (9) expresses as 

f, (l + £) m m! < "f 6 (l + f) m m! ~ 2"m! 

4^ (In 2) m m + 2 ~ ^ (In 2) m m + 2 ^ (In 2) m m + 2 ' { ' 

m—l m—1 m=n 1 / 6 +l 

When n (and so r) tends to +oo, the first sum on the right hand side of Eq. (10) 

converges to — — r— : — = 1.6702..., as a consequence of Lebesgue's dominated 
(In 2) m m+/ 
m>i v ; 

convergence theorem. And the second term goes to 0, as derived from Stirling formula. 
Finally, combining all these facts we obtain the announced lower bound on sj. □ 

The following Lemma 2 provides a lower bound on liminfpj. 
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Lemma 2 Let j = j(n) such that j(n) — > +00 anc? n/2 aJ — > 0. VFe /iaue 

liminf p 7 - > = .14846 . . . 

n 

Proof: Since (1 — x) > exp (— 2x) for x € [0, 1/2], we get the following 

£ow (m+1) 

fc=l 



v-!/n / n-l\Y m+1J ^ 1 / n \™+i / n(m + l)\ 

k=l v v 7 7 fc=l v 7 

El / n \m+i / n {ra + ] 

2™+i V2*J 6XP V 2^" 

fc=i v 



(m+1)! C/ m Mog 2 (n)J / : 

+ 7 .0 7= ,i + O ~ ,-, rr + 



2 m+i ( m + l) m + 2 ln2 (m + l) 2 4 m+1 \ v2 a3 + m + 1 n2 m+1 

By assumption, n/2 aJ — ► 0. So, the latter expression comes from Lemma 4, where the 
term U m (z) is defined in Eq. (18). Now, summing on m in Eq. (4) we derive 

t . > (_!]_, }_\ , m! su Pl/ |^(log 2 n)| 

J - I n j 2 m m m+1 In 2 ^ m 2 4 m 

m=l m=l 

By Eq. (19), numerical evaluation of the sums (with Maple) give liminf tj > .82092 . . . 

Finally, by using both lower bounds on Sj and tj (in Lemma 1 and in the above proof, 
resp.) we obtain the desired lower bound on liminf pj\ p\ = .14846. . . □ 

3.2 Proof of Theorem 1 

To prove Theorem 1, we use the following Lemma 3 that allows us to control the mean 
of the quantities of interests, whenever having a bound on the probability of success in 
each round. 

Lemma 3 Let (Xj)j>i and (i 7 i)i>i be two sequences of independent Bernoulli random 
variables, denoted by B(P,j) and B(Qi), respectively, and such that Pi > Qi for any i. 
By definition, 

P(Xi = 1) = 1 - F(Xi = 0) = Pi and P(Y; = 1) = 1 - P(y< = 0) = Qi. 

Let H = inf{j, Xj = 1} and K = inf{j, Yj = 1}, which may be regarded as a first success 
in each sequence Xi and Yi (resp.). Then, the "stochastic inequality" H <s K holds: 
for any non-negative integer k, P(H < k) > P(K < k). 
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Moreover, for any non- decreasing function f, 

E(f(H))<E(f(K)). (11) 

Proof: The first part of the Lemma can be proved by constructing a probability 
space ft in which the sequences of r.v. (Xi) and (Y) "live" and in which, for each 
to £ Q, Xi(uj) > Yi(uj). Since for any u, K(uj) > H(uj), the stochastic order is a simple 
consequence of this almost sure order on Q. Next, for any nondecreasing function /, 
f(K(uj)) > f(H(uj)) also holds almost surely and whence Eq. (11). □ 

We come back to the proof of Theorem 1. It is straightforward to check that, when 
n — ► oo, n/2 ab +1) — > 0. According to Lemma 2, if n is large enough, 

Pj > PiW+i ( 12 ) 

As a consequence, let N\ denote the number of rounds in Algorithm 1 and N{ = j*+G 
where G is a geometric random variable with parameter p*, then 

Nx < s N[, (13) 

(iVi is smaller with respect to the stochastic order than N[). 

We recall that the geometrical distribution with parameter p is given by p(l — p) k for 
any k > 1; this is the law of the first success in a sequence of independent Bernoulli 
random variables with parameter p. 

To show Eq. (13), notice first that Ni has the same distribution as inf{i, Xi = 1}, 
where the (Xi) are independent and Xi is -B(pj)-distributed. 

Next, N{ = inf{i, Yj = 1}, where Yj is ^(gj)-distributed for % = p\lj>j* +1 . Indeed, the 
first j* trials fail and afterwards, each trial results in a success with probability p*. 
Finally, Eq. (12) and Lemma 11 allow to conclude, 

E(N ± ) < E(N[) = f + 1/p* = log a log 2 (n) + 0(1). 

Let T\ = T\(n) be the time needed to elect a leader in Algorithm 1. Since N[ is 
larger than N\ for the stochastic order and since r ^ Yli=i \ al ~\ i s non-decreasing, by 
Lemma 3, 

/ jvi \ / K \ +00 j*+k 

e(Tx) = e < E D^l < E E( 1+tt, ')rf( 1 -rf)* -1 

\j=l / \j=l / fe=l i=l 

< C(pJ, a) log 2 (n) + O (log log n). (14) 

During a round, the mean number of awake time slots for a given station is at most 2. 
(Two at the end, one from the contribution of l/2 k , l/2 k being the probability to 
be awake at time slot k.) 

Since the number of rounds is smaller than, say: log a log 2 (n) + logloglog(n) with 
probability going to 1, we get the announced result. □ 
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Remark 2 It is easily seen that the algorithm and the convergence of the double sum in 
Eq. (14) (resp.) require the conditions a > 1 and a (I— p\) < 1, withp* = .14846 (resp.). 
The value of a may thus be chosen in the range (1, 1.743 . . .), so as to achieve a tradeoff 
between the average execution time of the algorithm and the global awake time. Thus, 
the minimum value of the constant C(p\, a) is C(p*, a) ~ 29.058 . . ., with a = 1.0767 

4 Analysis of Algorithm 2 

Two awake stations are needed in Algorithm 2: the one is only sending (the initiator) 
and the other is listening. The probability p'j that one station is elected in round j 
expresses along the same lines as in Eq. (4) with the corresponding probability. The 
probability of having exactly one initiator and one witness for k = j in a round is 

Hence p'j = t'j s'j , where 

^ = 2 9(*,n) i 1 
k=l 

All computations are quite similar to those in the proof of Theorem 1. Again, the value 
of p'j in Eq. (15) is obtained by asymptotic approximation (see Lemma 5 in Appendix A). 
First, 




Using this latter lower bound in Lemma 5 yields 

(°° -1 2m \ 

x - x - 1 n zm / nm\ \ 
m=l i=l J 

( _ ^ 1 (2m- 1)! ^ / 1 
> ex P o(l) - 2m +i 1 o — I — 3 sup K log 2 n 

\ m=l m=l s 

with the following value of the lower bound on s'- (computed with Maple): 

liminfs' > .19895... 

n J 



1(H^_1 

2 4? V Z> 



n-2 



Q(k,n)) 



and 



1=1 



Q(i,n))- 



(15) 
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Next, 1 — x 2 > e x when x is close to and so ( ] > — exp ( — ] ; whence i'- also 

\2j 2 V y/nj 3 

is bounded from below as follows, 



n\ n 2 ( 1 



fc=l m=l \ v 




Now, 1 — x > exp {—x — x 2 ) when x € [0, 1/2], and 

m — l I 1 V / 



m=l k=l 



n 


n 


2^ 


~ 4fc 


n 


n 


2^ 


~ 4^ 



oo I" I 1 / 2 \ m 

m=l fc=r 

where the last summation is starting from r = r(n) = 3/21og 2 (n). For such a choice of 
r, we have 2 r <C n <C 4 r (which is used in the following). 

Therefore, we can now use Lemma 5 to deal with the inner sum in Eq. (16) 



1 \ - I m n\ v - I n / n \ \ 

m=l \ v / £,_ r v 7 

n l/4 

e ^ 4r L <f iii „ 9 sup|K(log 2 n)| + o(l). 



m=l 

By Lebesgue's monotone convergence theorem, this sum converges to 

+oo 



( {2m-l)\ 1 \ 

m=l 



and the numerical value obtained is lim inf t'- > .39856 . . . 

Finally, from the lower bounds values on s 1 - and t'j, we find 

lim inf p' > .079294... 

n J 

Remark 3 The lower bound p 2 in Algorithm 2 is already defined in (2). Now, since 
p\ = .07929, a meets the new condition if it belongs to (1, 1.086 . . .), and the minimum 
value of the constant C{p^i «) is C(p 2 , a) ~ 52.516, with a = 1.0404 . . . 

Note also that Algorithms 1 and 2 can be improved by starting from k = fco, &o > 1 
(instead of 1) in the third line of the algorithms. Though the running time of each 
algorithm remains asymptotically the same, starting from k = ko > 1 reduces the awake 
time to (1 + e) log a log 2 (n) time slots for Algorithm 1, and (1.5 + e) log a log 2 (n) time 
slots for Algorithm 2 (with e = l/2 fc()_1 ). Yet, this also makes the running time longer 
for small values of n; whence the (obvious) fact that the knowledge of any lower bound 
on n greatly helps. 
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A Appendix 

The following two technical Lemmas are at the basis of the asymptotic complexity analy- 
sis of Algorithm 1 or Algorithm 2, in the proof of Theorem 1 and Theorem 2, respectively. 
They both use Mellin transform techniques [5, 6, 10] 

Lemma 4 Let r = r(n) such that, when n — > +oo we have simultaneously r — > oo and 
n/2 r — > 0. T/ien, /or a// positive integer m, 

v — *v / Tl\ m / TlTfl\ TX)} i / \ 

L(,2*J 'M"^) = m-+Hn2 + ^ (. '° 82 '"') 



fc=i 



Denote xe = 2i^7r/ln2. For any positive integer m, U m (\og 2 (nfj is defined as 

^m(^) = mm _! ln2 X] r ( m + »)exp [ - x<?ln(m)J exp(-2i&rz). (18) 

The Fourier series U m (z) has mean value and the amplitude of its coefficients does not 
exceed .024234. 

Proof: The asymptotic approximation of the finite sum in Eq. (17) is obtained by 
direct use of the properties of the Mellin transform [5, 6]. 
The sum 



+0O 

■ m 



k=i 



is to be analyzed as x — > oo. Its Mellin transform, with fundamental strip (—1,0), is 

T(s + m)2 s 



/*(«) 



m s+m (l - 2 s ) 



There is a simple pole at s = 0, but also simple imaginary poles at s = xe = 2^7r/ln2 
(for all non zero integers £), which are expected to introduce periodic fluctuations. The 
singular expansion of f*(s) in ( — 1/2,2) is 



/*(') 



T(m) I 
m m ln2 s 



— y 



r(m + x*) exp ( - ln(m)) 



m m In 2 ^— ' s — Y£ 



x \ m / mx 



\2k) ex P \ — ^F) -^( X )' one fi n< ^ s the results stated 
fc=i 

in Eqs. (17)) and (18). In Eq. (17), the periodic fluctuations are occurring under the 
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form of the Fourier series U m ^\og 2 (n)j with mean value 0. Besides, for any positive 
integer m, the maximum of the amplitude of the Fourier series is taken at m = 11 and 
it is rather small: 

Vm > \U m (z)\ < Y 2 m |r(m + y)| mA2M ^ 
1 Wl - ^ m m - 1 ln2 v ; 

eez\{o} 

The Fourier coefficients also decrease very fast (see e.g., [6] or [10, p. 131]). The error 
terms 0(2 m /n m ) and 0(n m /2 rm ) in Eq. (17) result from the "truncated" summation: 
1 < k < r. (The Mellin transform itself results in a O (n -1 ) error term.) □ 

Lemma 5 Again, let r = r(n) such that, when n — > +oo ; r — > +oo and n/2 r — ► 0. 
Then, /or aZZ positive integer m, 

El ( n 2 \ m / nm\ (2m -1)! 1 T ^ /, . .\ 

4^J eXp (-^J = 4^+1 In 2 + ^^l l0g2(n) J 
fe=i v 7 

+ (^) +0 (^)' <20) 

Again denote xi = 2i£7r/ln2. For any positive integer m, the above Fourier series 
V m (z) = - ^ m ^ m ln ^ F ( 2m + exp ( ~ 2i&rlog 2 (m)) exp (-2i£vrz) 

^GZ\{0} 

/las mean i>aZwe and the amplitude of the coefficients ofV m (z) cannot exceed 9.0054 10~ 5 . 

Proof: The asymptotic approximation of the finite sum in Eq. (20) is provided along 
the same lines as in Lemma 4. Now, 

. . x - 1 / x \ 2m ( mx \ 

9( x ) = 2^^{^k) exp V ~^J 

k=l 

T(s + 2m)2 s 



is to be analyzed when x — ► oo. Similarly, its Mellin transform is g* (s) - 

m s+Zm {l — 2 s ) 

with fundamental strip (—1,0). Again there is a simple pole at s = and simple 
imaginary poles at s = Xi = 2rf7r/ln2 (for all non zero integers £), which are also 
expected to introduce periodic fluctuations. 

As in Lemma 4 the singular expansion of g* (s) provides an asymptotic approximation 
of the finite sum in Eq. (20). The periodic fluctuations occur under the form of the 

Fourier series V m ^ log 2 (n)j . As U m (z), V m (z) has mean value and its coefficients have 
a very tiny amplitude. Their minimum is taken at m = 2, and 

Vm>0 \V m (z)\< y J r(2 ? + 1 ^ ) ' < 9.0054 10- 5 . 

4m m 2m-l l n2 

eez\{o} 

The Fourier coefficients also decrease very fast (see e.g., [6] or [10, p. 131]). Last, 
the "truncated" summation (1 < k < r) results again in error terms 0(2 m /n m ) and 
0{n m /2 rm ) in Eq. (20). □ 
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